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Abstract. We rigorously construct radial H 1 solutions to the 3d cubic focusing 
NLS equation idtip + Aijj + 2\ip\ 2 ip = that blow-up along a contracting sphere. 
With blow-up time set to t = 0, the solutions concentrate on a sphere at radius 
~ t 1 / 3, but focus towards this sphere at the faster rate ~ t 2 l z . Such dynamics were 
originally proposed heuristically by Degtyarev-Zakharov-Rudakov [2] in 1975 and 
independently later in Holmer-Roudenko [5] in 2007, where it was demonstrated 
to be consistent with all conservation laws of this equation. In the latter paper, 
it was proposed as a solution that would yield divergence of the L\ norm within 
the "wide" radius ~ || Vw(i)||7 j2 1/ ' 2 but not within the "tight" radius ~ ||Vu(t)||^2 > 
the second being the rate of contraction of self-similar blow-up solutions observed 
numerically and described in detail in Sulem-Sulem [131 Chapter 7]. 



1. Introduction 
Consider the 3d cubic focusing NLS equation on IR 3 : 

(1.1) id t i/) + A^ + 2|VfV> = 0, 

where ip = ip(x,t) £ C and x £ IR 3 , t £ R. The initial-value problem posed with 
initial-data in H.\ is locally well-posed and there is a unique solution in C([0,T); H£) 
on a maximal forward life-span [0,T). If T < oo, then lim^T || Vip(t)\\L2 = oo. 

During their lifespan, the solutions ip to (1.1) satisfy mass, energy and momentum 
conservation laws. 

(1.2) M(^(t))= [ \^(x,t)\ 2 dx= [ \Mx)\ 2 dx = M(^ ), 



(1.3) E(^(t)) = J (\V^(x,t)\ 2 -\ij(x,t)\ 4 )dx = E(ij ), 

(1.4) P(4>(t)) = J (^V^ - dx = P(if> ). 

The case T < oo occurs for a large class of initial data (see the discussion in Holmer- 
Platte-Roudenko [6]) and in this case we say that the solution blows-up in finite time. 
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Numerical results (see Sulem-Sulem [131 Chapter 7]) describe the existence of self- 
similar radial blow-up solutions of the form 



(1.5) 



u(x, t) 



W) u (m) ei ^ T ~' )w '' thKt)=vm ^ T) ' 



where U = U(x) is a stationary profile satisfying the nonlinear elliptic equation 



1.6) 



AU-U + ib{U + y ■ W) + 2\U\ 2 U = 0. 



Numerics suggest that b > can be selected so that a nontrivial zero- energy solution 



U to (1.6) exists. Asymptotics as |x| — > oo show that any solution of (1.6) fails 



to belong to Hx (and hence H].). Nevertheless, it is expected that H\ blow-up 



solutions to (1.1) well-modeled by (1.5) near the origin exist, and the proof remains 
an important open problem. 

In this paper we consider only radial solutions, thus, P(tp) = 0. We construct a 



family of finite-time blow-up solutions to (1.1) with different dynamics - they focus 



toward a sphere, while at the same time the radius of sphere is shrinking. We study 
the problem by converting it (via the time-reversal, u(t) i— > u(t), and time-translation 



symmetry of (1.1)) to one with solutions defined on (0,to] f° r > that "start" at 
to > and blow-up at t = as they evolve backward in time. 
Our main result is the following: 

Theorem 1.1. For any e G R there exists a radial solution ip G C((0, to], iT 1 (lR 3 )) 



of (1.1) with E(ijj) = e, which starts at time to > and, while evolving backwards in 
time, blows-up at t = 0, and has the following form: 

(1.7) ip(r,t) = e im e iv{t)r/2 X{t) <p(\{t){r - q{t))) +h(r,t), r = \x\,x G M 3 , 
where <p(r) = sech(r) and the time- dependent parameters (q, v, A, 9) satisfy 



(1.8) 

and 
(1.9) 



lit) 


= qot 1/3 


% 


_ 2V331/6 


X(t) 


= \ot- 2/3 


A 


- 2 _2 / 3 3 -1 / 3 


v(t) 


= vor 2 " 


v 


— 2 1 / 3 3~ 5 / 6 


9(t) 


= e r^ 


Oo 


- 2~ 1 / 3 3~ 2 / 3 



\\h(t)\\ LHm +t 2 / 3 \\h(t)\\ mm +r 1/3 |||a;| h(t)\\ L2m < Ct 1 / 3 . 



Note that it follows from (1.8) and (1.9) that 

\\Vip(t)\\ L 2 m ~r 2/3 and \\\x\^{t)\\ L 2 {m ~t 
Remarks. 1. Let <p be a solution of the nonlinear elliptic equation 



1/3 



;i.io) 



(p + Aip + 2\(p\ 2 (p = 0. 
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The unique (up to translation) minimal mass H 1 solution of (1.10) is called the 



ground-state. It is smooth, radial, real-valued and positive, and exponentially decay- 
ing (see, for example, Tao [HI Apx. B]). Note that tp(x,t) = e lt (p(x) is a solution of 
( 1.1[ ). In one dimension, i.e., d — 1, the ground state of ( |1.10[ ) is explicit, namely, 
(p(x) = sech(x). It is this function we use for the profile of the blow-up in (1.7) (and 



not the 3d ground state solution of (1.10)). 



2. It follows from (j 1 . Sp that the solutions given by Theorem 

MU 2 (iR 



1.1 



satisfy 



L 2 



= \/2. Applying scaling, one then obtains contracting sphere blow up solu- 
tions with an arbitrary mass. 

3. We do not address in this note the questions of stability of these solutions (the 



numerical and heuristical results of [3] indicate that the behavior (1.7), (1.8) should be 
stable with respect to radial perturbations). Still, let us mention that the arguments 
developed in the proof of Theorem can be easily adjusted in order to produce 
blow up solutions of the form 

ii>(t) = S(t) + f + o m (l), *->0, 

where S(t) is a solution of ( 1.1[ ) given by Theorem |1.1| and / is an arbitrarily smooth 
radial decaying function vanishing to a sufficiently high order at the origin, very much 
in the spirit of [I]. 

4. A numerical visualization of the blow-up described by Theorem LI is given in 
Figure 1, courtesy of R. Platte, Arizona State University. 





FIGURE 1. A contracting sphere collapse in the 3d cubic NLS: (a) a con- 
centrating layer of the collapsing solution is around r ~ 0.6, (b) the eventual 
collapse at the origin. Courtesy of R. Platte, Arizona State University. 



The paper is organized as follows: In £|2j we construct (Prop. 2.1) for arbitrary 
6 N, an approximate solution ip^ N ' satisfying 



(1.11) 



(N) 



4 

where 



JUSTIN HOLMER, GALINA PERELMAN, AND SVETLANA ROUDENKO 



t 2k/z \\V k H 



(N) 



L 2 (R3) 



< fc t (27V - 1)/3 , 0<t<t (N), 



with implicit constant independent of N. To accomplish this, we set 
(1.12) V (A °M) = e ie e ivr/2 \U {N \p,t), 

where p = X(r — q) is the recentered and rescaled spatial variable and q, v, 8, A are 
t-dependent parameters. We then seek to suitably define the parameters and profile 
C/W. We impose the a priori conditions 



Q 



which achieve convenient simplifications of the equation that arises when (1.12) is 
substituted into (1.11), and enable us to solve for v and 9, once q and A have been 
>r add 

U( N > — cp (recall (p(p) = sechp). Then (1.11) will follow from 



specified. For additional convenience, we set u> = A 1 q 2 (effectively replacing A by 
u) ancQx 



;i.l3) mHx + £(q,oj)x + F(q,u) + F^(jl))\<Ce >fl t 



f (27V+3)/3 



< jJL < 1, 



where TL is the ID cubic linearized Hamiltonian - %, £, J 7 , and F- 2 are defined 



precisely in g2| We then seek a solution to (|1.13|) by expanding 
(1.14) 



N 



N 



27V+2 



k=0 



k=0 



k=l 



Substituting the expansions (1.14) into (1.13), we can organize terms by their order 
in t and carry out a recursive construction. The parameter coefficients and Uk are 
selected at each stage to satisfy solvability conditions for Xk+i- 

In ^3j we consider the exact solution ip with data ip\t=e — ip^(s) for any < e < 
t (N), and prove (Prop. 3.5) that for N sufficiently large, the comparison estimate 

2 

(1.15) £ t 2k ^\\ V\m - ^ N \t))\\ LHm < t 2N " , e < t < t (N), 

k=0 



and localization bound (Prop. 3.6) 

(1-16) ||^(t)|| L2(R 3) <t 1/3 



hold, independently of < e < t (N). The estimates (1.15) and (1.16) are proved 
using a Lyapunov functional G(h) (see (3.8)), built from energy, mass, and radial 
momentum. 



Finally, in §4, we complete the proof of Theorem 1.1 using estimates (1.15), (1.16) 
and a compactness argument. 



1 In fact, once x is constructed, we take to be a cutoff of ip + %■ 
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This result was announced and described in detail 0| by G. Perelman in the Anal- 
ysis seminar at Universite de Cergy-Pontoise in December 2011. In February 2012, 
Merle, Raphael, & Szeftel posted a paper [Sj on arxiv.org proving a closely related 
result. We stress that our result was obtained independently of [8]. 

1.1. Acknowledgements. J.H. was supported in part by NSF grant DMS-0901582 
and a Sloan Research Fellowship (BR-4919). S.R. was supported in part by NSF 
grant DMS-1103274. This work was also partially supported by a grant from the 
Simons Foundation # 209839 to S.R. 

2. Construction of an approximate solution 



We start with writing a solution to (1.1) as 



(2.1) ip(r,t) = e im e iv{t)r/2 \(t)U(\(t)(r-q(t)),t). 



Substituting it into (1.1) and introducing a rescaled coordinate 

p = \(r- q), 

we obtain that U = U(p, t) solves 

2 U p -i^(l + pd p )U- i v u. 



p + Xq P A 3V r P! X(p + Xq) 

Imposing the following two conditions 

v(t) = q'(t), 

(2 ' 3) m + ^ + ^M = x(tf , 



we obtain 



1 ,„,,„ 2 



(2.4) i — U t + U PP -U + 2 \U\ 2 U = — U, 

X z p + Xq 



We next denote 

(2.5) u 



Xq 2 
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and rewrite (2.4) as 

(2.6) iqWUt + Upp-U + 2\U\ 2 U 



■ 2qu} U p 

pcoq + 1 



+ rfq^p U + t{2q'q 3 u 2 + u'uq 4 ) (l + p 8 P ) U 



qq 3 u 2 

', 

puq + 1 



U. 



Define the Pauli matrices by 



<7\ 



1 

1 



and (T3 



1 
-1 



Split U into the localized part cp(p) = sechp and the external part x(p^) : 
(2.7) U(p,t) = <p(p)+ X (p,t). 



Substituting (2.7) into (2.6) and using the vector notation \ 



we obtain the 



following equation for x' 

(2.8) Hx + l(q, u)x + F(q, u) + F^\x) = 0, 

where the components of the left-hand side are defined as follows: 

• the Hamiltonian corresponding to the linearization of the ID cubic NLS idtu + 
d 2 u + 2\u\ 2 u = around the soliton u(t, p) = e lt sech p: 

U = (-d 2 p + 1) (7 3 - V(p)a 3 - 2<p 2 (p)a 3 a l , 

• the linear operator 
A. ,2a 2 <?W 



l(q, oj) = —iq oo d t 



pooq + 1 



d p a 3 + i(2qq 3 uj l + u'ujq 4 )(l + pd, 



qq 3 uj 2 

', 

puq + 1 



• the terms independent of x~- 

F{q,uj) = l(q,u)<p 

and lastly, 

• the higher order terms grouped into 



F^(x) = -4^|x| 



1 

-1 



2ip 



X 

-x 2 



21x1 



X 

-x 



Proposition 2.1. There exist coefficients (<jJk)k>o, {<lk)k>o> with qi arbitrary and 
u = \, and smooth exponentially decaying functions (Xk{p))k>i> 



(2.9) 



^ p[ d l p Xk(p)\ <Q >tM;k , pe 
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for any k > 1, I > 0, and < p < 1, verifying the orthogonality conditions 
(2.10) (Rex/cV^L^R) = (Imxfc,y?)L 2 (R) = 0, V7c > 1, 

such that the following holds. For any N > 2, if we define 

N N 

u(t) = ^ N \t) = 1 + $>,i 2fc / 3 , q(t) = q {N) (t) d ^Y,<lkt i2k+1)/3 , 

k=l k=0 
2N+2 

7 / , \ rm i , \ dcf 

anc 



k=l 



then (2.8) approximately holds: 

(2.11) \d l p (Ux + l(q, u)x + T(q, u) + F^ 2 {x)) \ < Qj^^e'^, 

for < p < 1, p > -\q G l t~ l/z , < t < t (N), with some t (N) > and with C t ^ 
independent of N. 

N N 2N+2 

Proof Writing u(t) = J> fc t 2fc / 3 , q(t) = ^ g fc t (2fc+1)/3 , X (p,t) = £ t k '\ k {p), 

k=0 k=0 k=l 

uq = 1 and substituting this ansatz into the expression 1-ix + ^(<L w )x + ^(q^) + 
F^ 2 (x), we get 

2N+2 

Ux + l{q, uj)x + Hi, w) + F^{x) = Yl tk '^ U Xk - Ac) + S(t), 

k=l 

S(t) = ~ Yl t kri D k + h{q^)^{\)+X), 

k>2N+3 ^ ' 



(2-12) 



where 



N 2qcj(qujp) 2N+ \ q'q 3 uj 2 (qtjp) 2N+3 

h(q,u) = n d p a 3 + i — . 

puq + 1 pooq + 1 



Df \ being contributions of F{q,u) —h(q,u)(p ^ j J, (l(q,u) -h(q,uj))x 

and -F- 2 (x) respectively: 



fe>i 



(/(g,a;)-/ 1 (g ) a;))x = -^t fc/3 J Di 1) (p) ) 

fc>2 

^- 2 (x) = -E^ 2) ^)' 



fc>2 
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D%> = ( * w ). Note that all these sums contain only finite number of terms. 

The following structural properties of will be important for our analysis. 

(i) Dependence on qj, 00 j xf depends on qj, 00 j, j < (k — l)/2, and i = 1,2, 
depend on qj, Uj, j < (k — 2)/2, and on \pi P < k — 1, only: 

Df = D^(p ] q v u J ,j<(k-l)/2), 

Df = Df{p- qj,uj,j <(k- 2)/2; X P ,P<k-l), 1 = 1, 2. 
In addition, one has 

(2.13) t^D^ = 2q (t)<p P - tq' (t)q 3 (t)(if + 2p<p p ) - ^(t)q%(t)pp, q (t) = q t^ 3 , 
and for k = 21 + 1, Z > 1: 

(2-14) flS^^ + cS, 

where -D^+i = ^a+iG ; Qji^jiJ — I ~ 1) depends on qj, Uj, j < I — 1, only and -D^+i 
is given by 

(2.15) 

t (2Z+l)/3£)(0,0) _ 

2(ffi(*) + g (*M(OK - ^ WW + 6g '(i)g 5 (*k(*) + 3g '(t)g 6 (tM(t) W 
- % 3 (M*) + 3g (£)g 2 (t)g,(t) + 2g (t)g 3 (£)u^))(<P + 2p<p p ) - zg 4 (f)u;;(t)(<p + p<p p ), 

«(t) = ^ (2W) / 3 , ^(t)=^ 3 . 

(ii) Smoothness and decay: if all Xji j < k — 1, verify (2.9), then a similar estimate 
holds for D k : 

(2.16) \e M d l p D k (p)\<Q^ kl pGM, 
for any I > and p < 1. 

(iii) Parity: if for any j < k — 1, Vj d = KeXj and Uj = f ImXj satisfy 

(2.17) Vj (-p) = (-l)^(p), %(-p) = (-l) i+ S(p), 
then one has the same property for Q£ = f Re/^, ^ = f ImDk : 

(2.18) #(-p) = (-i) fc £ fc + (-p), g k (- P ) = (-i) k+1 g k (- P ). 

Our goal now is to show that we can find coefficients go, Qi, ■ ■ ■ , gjv, ^i, • • • , wjv (go 
being as in Theorem |1.1| and q2 arbitrary) and functions xu ■ ■ ■ > X2N+2 satisfying (2.9), 
plQj ), J2T7| ) such that 

(2.19) nxk = D k , k = l,...,2N + 2. 
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Estimate (2.11) will then follow from (2.12). As we see below, the equations with 
k = 5 require some special care, therefore, we will go explicitly through k — 1, ... ,5 
and then we will proceed by induction. 

It will be convenient for us to rewrite (2.19) as 

(2.20) L+v k = g+, L-u k = Gk, k = l,...,2N + 2, 

where L + = —d 2 + 1 — Qip 2 , L_ = —d 2 + 1 — 2ip 2 are self-adjoint operators in L 2 (IR) 
with the spectrum cr(L+) = {—3, 0} U [1, oo), cr(L_) = {0} U [1, oo). One has 

L_<p = 0, L + ip p = 0. 



For k = 1 the equations (2.20|) give 

(2.21) 
with 



L + vi = Gt, L-Ux = Q x 



see (2.13). 



Since (Q x ,tp) L *{R) = + 2 PW V?) 



0, the equation L_u x = Q x has a 



unique solution u x in L satisfying (iti, <£>)l2( K ) = 0. One can compute it explicitly: 



o 



\f*p\\h( 



\<p\\h( 



The solvability of the equation L + v x = Q x is subjected to the condition 

{Gi, <^p)i 2 (K) = o, 

which gives 

go - Aq 7 Q = 0, A 

Therefore, we find 

(2.22) g = 12 1/6 . 



W\\ 2 m 



36||^|| 



1 

12' 



Under this solvability condition, the equation L + v\ = Q x has a unique solution v\ 
satisfying (v 1: (P p )l*{r) = 0, which is given by 

{d p {p 2 <f),<£ p ) L 2 m 



Vl 



go 



(d P (p <p) + cxtpp), ci 



||^p|U 2 (M) 



Note that v% is an odd and u% is an even function of p, in agreement with (2.17), and 
they both satisfy (2.9). 

Consider fl2~2u) with k = 2. We have 

(2.23) L + v 2 = Gt, L_u 2 = G 2 i 



10 JUSTIN HOLMER, GALINA PERELMAN, AND SVETLANA ROUDENKO 

where Gt = Gt(p\ %)■> $2 = G^iP'i 9o) are smooth exponentially decaying functions of 
p, verifying (2.16), (2.18). Since Gt is even and Q% is odd, the solvability conditions 

(^2" 5 y)i»(M 



are satisfied and (2.23) has a unique solution t> 2 , w 2 



with t> 2 even and it 2 odd, verifying (2.9). 
Next consider k = 3,4. For k = 3 we have 

(2.24) L+w 3 = <#(go>?i,Wi,Xi,X2), L_u 3 = G^(q ,qi,u 1 ,Xi,X2) 



with g ^(g , g i, xi, x 2 ), £3 (<?o,gi,wi,xi,x 2 ) verifying ( |2.16[ ), ( |2.18[ ). Using ( |2.14[ ), 
(2.15), (2.22), one can write Gt , G3 as 



(2.25) 



Gt 



2 (ft. + q uJi)<f p + 4(2gi + g ^i)p<P + 0; 



3 • 



2 2 



where ^ = Gt (<?o, Xi, X2), G 3 = G 3 (qo,XuXi) depend on q ,Xi,X2 only, and there- 
fore, are determined by now. 

Using (2.25), one can write the solvability conditions {G% , (p)l 2 (m.) — (Gt > <Pp)l 2 (w£) — 
as 



qt\W\\h< 



-Wi + (C/ 3 , V5)l2( R) = 0, 



2{q x + g ^i 
which gives, 



2 

l 2 ( 



2 (2gi + g wi)l|<p|li2 (M) + (£3^) 



to 1 



3 (^3 i^L^ 
9ollvlli2(K 



L 2 ( 



qi 



-5 qoUl+ 10 



(£3 > <Pp)-L 2 (R) 



2 



With this choice of U\, qi, (2.24) has a unique solution t> 3 , u 3 verifying (2.9), (2.10) 
with t> 3 , m 3 being of the same parity as Gt, V3 is odd and U3 is even. 
The case k = 4 is similar to that of k = 2: we have 



(2.26) 



£+i>4 = Gt, L_Ui = G, 



4 ) 



with £4 even and C/J odd, depending on q , q ly ui, Xi, X2, X3 only. Thus, (2.26) has a 
unique solution t> 4 , w 4 with t> 4 even and w 4 odd, verifying (2.9). 
Consider the case k = 5. We have 



(2.27) 



L + v 5 = Gt, L_u 5 = G\ 



5 , 



with^ = ^(g ,gi,g2,wi,W2,Xi, •• • ,X±) odd and G 5 = G 5 (<?o, ft, ft, wi, w 2 , Xi, • • • , X4) 



even. 
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It follows from (]2T5]), ( p^2[ ) that 
(2.28) 



£ 5 + = 2(q 2 + q uJ2)Pp + g(?2 + 3g w 2 )p<p + £5", 



2 4 
05 = -gt 4 ^ + go w 2)(<P + 2p<p p ) - g<?oW<^ + P<Pp) + Gs, 

where Qf = G^(qo, <Zi> ^i>Xi> • • • , Xa) are determined by now. Substituting this rep- 
resentation into the solvability conditions (G5, <p)l 2 (r) = 0, (^ 5 h , v 9 p)l 2 (ir) = 0, we 
obtain 

( ^ 2 ~ 2a 4 IMI 2 ' 

zc /0ll ( PllL 2 (IR) 

(2-30) - 4q u 2 y p \\ 2 L^R) + <P P )vp 

which means that one has to have 
(2-31) QoiGt, ¥?p)l2 (M ) = 2(£ 5 ~ (f) L 2 { 



0. 



Note that (2.29), (2.30) do not contain q 2 , therefore, it can be chosen arbitrarily 
(as soon as we can show that (2.31) holds). To check (2.31), we proceed as follows. 
Consider = 1 + u.t 2 / 3 , 9 «(t) = got 1 / 3 + Ql t, X {1) (p,t) = £j =1 t k ' 3 X k(p)- By 

our construction, they verify 
(2.32) 



1 + w (i) g (i)p 



rr«s(p,f), 



where = (_^), = —^5" — ^5 , and S(p,t) admits the estimate 

(2.33) \e^d l p S{p,t)\ <QJ 2 , peR, 
for any I > and p < 1. Define 

(2.34) z(x,t) d = f e i6 ^ +iv ^ 2 X 1 (t)(^ + x (1) )(Ai(t)(r - (*)),*), 

where Ai, Vi, 61 are given by (2.3), (2.5) with q = q^ and to = co^. Then z{t) solves 

(2.35) iz t + Az + 2\z\ 2 z = Ro{x, t) + Ri(x, t), x € M , 
where (by fl2T32| and fl2T33~] )) 

i2o(x,t) = -t^e^^+^^A^O^Ai^^-g^Ht))), 



and i?i satisfies 
(2.36) 



I-Rlll i 2 (R3) < C t 2 ^ 3 . 



Consider the energy E(z(t)) = J(\Vz(x,t)\ 2 — \z(x, t)\ 4 )dx. It follows from (2.35), 

d 



(2.36), (2.34) that 



(2.37) 



dt 



E{z{t)) = at- 1 + 0{t- 2 / 3 ), t^0, 
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with 
(2.38) 



JUSTIN HOLMER, GALINA PERELMAN, AND SVETLANA ROUDENKO 



a = — - 



<?o(^5 + >V?p)l2(M) - 2(£ 5 ,<f) L 2 { 



On the other hand, using the definition of z(t), we can write an expansion of E(z(t)) 
m powers of *V3 : 



(2.39) 



E(z(t))=Y,t k/3 c k + 0(t 1 / 3 ), t^O, 



fc=-4 



with some constants e&. Comparing this expansion to (2.37) one gets c„4 
c_i = and a = 0, which is precisely (2.31). 



Verifying the solvability conditions (2.29), (2.30) results in the fact that the system 



(2.27) has a unique solution v 5 , u 5 satisfying (2.9), (2.10) with v 5 being odd and w 5 
being even. 



To finish the proof of Proposition 2.1, we proceed by induction. Suppose we have 
solved (2.20) with k = 1, . . . , 21 — 1, I > 3, and have found qj,Uj, j < I — 1, and Xpi 
p<2l-l, verifying Q, fl2lo] ), and §2~17\ . Consider k = 21, 21 + 1. For k = 21 we 
have 



(2.40) 



L+v 2 i = QtMi,u h i <l-l;x P ,P<2l- 1), 
L-U21 = <3 2l (qj,Uj,j < I - l;x P ,P <2l- 1), 



with Q^i being even and Q^ x being odd, verifying (2.16). Therefore, the equation (2.40) 
has a unique solution V21, U21 with t> 2 ; even and U21 odd, that satisfies (2.9). 
Consider k — 21 + 1. We have 



(2.41) 



L_u2i+i = Gu+iiVj^jJ < l;x P ,P < 2i), 



where Q21+1 * s °dd and {? 2 ;+i ^ s even - From (2.15), (2.22) one has 

4 



(2.42) 



S21+1 = 2 (« + <?cmV p - g((2Z - / - 7)^ - 3q wi)pip + 

2g 3 2lq 4 
G21+1 = — ^-((l + 2) ft + goWj)(v + 2p^ p ) ^-%(<p + p<Pp) + Q 2 i + i, 



where Q 2 i+\ depends only on qj,ujj,j < I — \\\vV — Using this representation 
one can write the solvability conditions (Q 2 i+ii v)l 2 {r) = 0, (G21+H Vp)l 2 (r) = as 



q^lui 



\<P\\hcR) + 



0. 



2 qi (2l 2 - l-6)\\(pp\\ 2 L 2 {M) - 4qaWi\\<Pp\\h(&) + ($21+1, ^p)l 2 ( 
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which gives 



4 (^2«+1^)l2(K) (Gtl+V C Pp)L 2 ( 



2(2/ 2 — I — 6) y ^Zoll^plHanR) Ir/'ll/,^, 



Ml? 



After verifying the solvability conditions, one finds v 2 i+i, u 2 i+i as the unique solution 
of gjjt satisfying (gig. □ 



Define truncation functions G C£ 



(2.43) flo(C) 
For each iVeN, introduce 

(2.44) 

Note that 
(2.45) 



o, |CI> 1 ■ 



200 
1 

100 



by 
0i(O 



1 for |C| <h 
for|CI>^ 



(2iV+2 \ 



By construction each approximation solves (see equation (2.4)) 

1 „fj\n _.vaa „/nn _ ,../nn,o„/*A 2 



(2.46) 



i _L [/f ) + c/W _ [/W + 2 |[/W| 2 C/W 
A 2 w 



P + Ag p 



+ 



2A 3 



A' 



UW-i^(l + pd p ) uM-i 



q' 



X(p + \q) 



where the remainder 1Zn is supported in {\p\ < 10( ^ } by (2.43) and (2.45), by (2.11) 
for each N e N there exists a time to(-^0 > such that 7Zn satisfies the following 
bounds 

(2.47) \d k p U N {p,t)\ < C k e-" M t {2N+3)/3 for < t < t (N), k = 0, 1, 2, 
with the constants p and C\. independent of N. 



Observe also that from (2.44) (the first term is time- in dependent <p), we have by 



( |2T9[ ) that 

(2.48) \d k p U {N \p,t)\ < Ce-^ for < t < t {N), fc = 0,l,2, 
For each JVeNwe set ip^ by 

(2.49) ^ N \r,t) = e ie{t) e iv{t)r/2 \{t)U {N) {\{t)(r-q(t)),t). 
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Then solves 

(2.50) id t ^ N) + M {N) + ^ {N) \ 2 ^ {N) =iH {N \x,t), 
where 

(2.51) H^ N \r, t) = ie i0 Wj v ® r / 2 \ 3 K N (\(r - q(t)),t). 



By ( ]2.47[ ), it follows that 

11^^(011^3) < C fc A 3 +^< 2 ^+ 3 >/ 3 ||l |r -„ (7| < A -x]| z , 2(M 3 ) , 

and hence, 

(2.52) t 2k ^d k r H( N \t)\\ L2(R s ) <C k t( 2N - 1 ^ forO<t<t (iV), fc = 0,1,2, 



Proposition 2.2. Given e E M., we can choose q 2 in the construction in Prop. 2.1 
such that the energy of the approximate solution E(ip( N ') estimates as 



(2.53) 



|£(V> (JV) )(t) - e| < Ct^- 1 for0<t< t (N). 



Proof. The following two estimates are straightforward consequences of (2.50) and 



(2.52) 



(2.54) 

and 
(2.55) 



d 
dt 

d_ 
dt 



U {N \t)\\h m 



<Ct— for < t < t (N), 



<CX 2 t— for 0<t <t (N). 



Consider q- 2 (t)\- 3 (t)E(ij^)(t). By construction of V (A °, it admits an expansion m 
positive powers of t 1//3 (compare to (2.37), (2.38)): 



K(N) 



(2.56) 



q- 2 (t)\- 3 (t)E(^)(t) = ^t k/3 + 0(t°°), t 0, 



A:=0 



where (3^ depends on Xp, V < k, and on ujj,qj, j < k/2. More precisely, for k = 21 
one has 

4 



(2.57) 



#H = i(wi + (2Z + 3)^-) + Aa, Z>0, 



with /3 2 ; = $2i(Xp,P < 2l',Uj,qj,j < I — 1) depending on x P , P < 2/, and on u/j,g 3 -, 
j < Z — 1, only. From (2.55), (2.56), we deduce /3o = ■ ■ ■ = /?3 = and 

|£(^ (A °)(t) - <? 4 /3 4 | < Ct 2 ~T- 1 for < t < t (JV). 



Therefore, if we fix q 2 by requiring |(w 2 + 7^) + /?4 = eg^, then we get (2.53). □ 
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3. Comparison to exact solution 

Let us begin with some notational conventions for this section. We exclusively use 
the real inner product 



(3.1) (ipi, ^2)l2(r3) = Re y ^i(x)i) 2 (x) dx 

on L 2 (IR 3 ) as well as the related ID version 



( , 0i,V'2}x2(o<r<oo) = Re / ViO")^(r)dr. 

Jr=0 

Note that 

(^1,^2) V(b?) = 47r(r^i(r),r^ 2 (r))L2(o<r<oo)- 
For a (densely defined) functional A : L 2 (IR 3 ) — > R, via the inner product (3.1) we 
identify A'(if)) with a function and A"(tJ)) with an operator, which is self-adjoint with 
respect to (3.1 ). Let 



and note that, with respect to the inner product (3.1), J* = — J. Define the Poisson 
bracket 

{A,B}^) = (JA'(^),B'^)), 

which yields a functional {A, B} : L 2 (IR 3 ) -> E. 

We find it convenient to state estimates in terms of the time-dependent scaled 
Sobolev norms: 

k 

(3.2) \\h\\ xkm = J2 t2l/3 \\^ h \\^ 

Most frequently, we use the case k = 1: 

II^IU^k 3 ) = ^ 2/3 j| V/i|| i 2 (]R 3 ) + ||/i||l2 ( 

By default X 1 = X 1 (IR 3 ), although we shall have occasion to use the variant X x (0 < 
r < oo), which has the expected definition. 

We shall frequently need the following radial Sobolev inequality. For any radial 
function /, 

(3.3) ||/ • W/»IIl-C*) S ^llWIli^) ll/llK( K 3) < t-^\\f\\ xHM sy 

Recall the function ipW(t) defined on < t < t (N) by ( |2.49| ). For any < e < 
to(N), let be the solution to (1.1) such that 



(so ^ depends on e and iV but this is suppressed in our notation). Let 
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Each of the Propositions 3.2, 3.3, 3.4, 3.5, and 3.6 in this section references ip and 
h defined as above. All bounds stated will be valid for < t < t (N). All implicit 
constants (indicated through the notation <) will be independent of e and N. 
Note that (1.1) can be written as 

(3.4) d t i) = JEfty), 



while (2.50) can be written as 

(3.5) d t ^ = JE'{i,^) + H^ N \ 



where H^ N ' is defined in (2.51) and satisfies (2.52). 
Recalling (2.43), define the radial localized momentum 

2 



(3.6) V q ty) 
Define the functional 

(3.7) W{if>) = I 1 



Im 



v 

4A2 



0i 



r — q 

q 



ijj d r ip dx. 



O(l) 



y 2 W)+ v2 

o(t 2 / 3 ) 



1 

0(t 4 / 3 ) 



where M(tp), V q (ip), and E(ijj) are the mass, localized radial momentum, and energy 
functionals defined in ( 1.2[ ), (3.6), and ( 1.3[ ). Define the Lyapunov functional 

(3.8) G(h) = W{$) - W(^ N) ) - (W'(^ N) ), h). 

The statement and proof of the following Lemma are based on Holmer-Lin [U 
Lemma 5.1]. 



Lemma 3.1. For G{h) as defined in (3.8), we have 

(3.9) d t G(h) = {E, W}(i>) - {E, W}{^) - ({E, W}'(^), ft) - ^ + £ 2 , 
where 

E x d ^ f (W"{ipW)H( N \ h) + (W (</> {A °), J&ty) - J&tyW) _ jE"(ipW)h), 
£ 2 is the result of the time derivative landing on any of the parameters in W , and 

H (N) 

is as in ( |3.5[ ). 
Explicitly, 

(3.10) £ 2 = ( l -vv\- 2 - ^ 2 A- 3 A)(M(^) - M(^ (JV) ) - (M'(^ (A0 ), h)) 

+ (-i)\~ 2 - 2v\-*\){VM - W iN) ) - (P^ (N) ), h)) 
-X- 2 X(E(i;)-E(^)-(E'(^),h)). 
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Proof. We write expressions for the time derivatives of each term in (3.8), dropping 
the terms that lead to (3.10). First, from (3.4), 

(3 ' U) = (W'&),JE'&)). 



Next, from (3.5), 

d t w(tfj iN) ) = (w\^ N) ),d t ^ N) ) 

Finally, from ( |3.4 ) and (3.5), 

^(W(v (JV) ),/i) = (w"(^ N) )d t ^ N \h) + (w\^ N) ),d t ^-d t ^ N) ) 

(3.13) = (W"(^ N) )(JE'(^ N) ) + H {N) ), h) 

+ (W'(^ N) ), JE'(tfj) - JE'{^) - ffW). 



Taking (3.11) minus (3.12) minus ( [3~13| , and noting the cancelation of (W(^ (A °), # (A °) 
in (3.12) and ( 3.13[ ), we obtain 

(3.14) d t G(h) = (W'W, JE'^)) - (W'&W), JE'{^)) 

- (W"(^ N) )(JE'(^) + HW), h) 

- (W'(^ N) ), JE'(iP) - JE\^ N) )). 

Since 

{E,W}(iP) = (JE\i;),W'(4>)), 

we compute 

(3.15) ({E, W}'(^ N) ), h) = (JE"(^ N) )h, W'(ip w )) + (JE\^ N) ), W"{^ N) )h). 



Now, (3.9) follows from (3.14) and (3.15). 



□ 



Recall the definition of the X norm in (3.2). By processing the terms in Lemma 



3.1 we obtain the following upper bound on dtG: 



Proposition 3.2. Let G(h) be defined as in (3.8) and suppos^j H^Hx 1 < 1- Then 
(3.16) \dtG(h)\< ^111111^+ t\\H^\\ 2 xl . 

This bound is valid for e < t < to(N), and the implicit constants are independent of 
e and N . 



This is part of the bootstrap hypothesis. In fact, the argument ultimately yields the much 
stronger estimate H^Hx 1 < t 2N ^ 3 . 
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Proof. To apply Lemma [37L| we need to bound Si, S 2 , and to compute the Poisson 
bracket {E, W}. 

Computation yields (here, 81 is short for 9\{(r — q)/q)), 



(3.17) 



= -2d 2 tp - Ar^drip - 4|-0| 2 V>> 
and, for convenience, taking a = d r {9\) + 2r _1 # 2 , we have 



(3.18) 



First, we address E\. We have 



(3.19) 



2(1 + 



\v 2 \- 2 ) ^ - i v\- 2 {29%^ + aif) W) 



+ A-^(-2A^)-4|V (7V) | 2 V (7V) )- 
e>(t 4 / 3 ) 



By (3.19) and (3.3), we obtain 

\\w'(^ N) )\\ L2 <u w \\ x ?<i. 

However, this ignores the fact that the leading order term in is the push-forward 
of ip(p) = sechp, which leads to a vanishing term in (3.19). Thus, we in fact have 

(3.20) \\w\^)\y<t x i\ 

We have 

(3.21) E?W) - E'{^) - E"(i>W)h = tfj^h 2 + h 3 

(ignoring complex conjugates). In the pairing (W'(ip^), JE'(ip) — JE'(^ N ^) — 
JE"(ip( N ^)h), each term will have a factor of tp( N \ which localizes to r ~ t 1 / 3 . 
By fl3~21~j ), Q, and fl3~20] ), it follows that 

\{W'{^ N) ), JE'(iP) - JE'(^ N) ) - JE"(^ N) )h)\ 

(3.22) < ||M^'C^ CJN0 )IIj^CII^ CJ>0 II^~I1^ " 1, — *x/a|| J &~||/*||^ -H ||/«- - 1, — *^3|lx,~ll^llx^) 

From d3T9"l), we find that W"(ip^) is an operator of the form 

(3.23) W"(^ N) ) = 0(1) + 0(t 2/3 )(d r + a) + £>(t 4/3 )(A + (^ N) ) 2 ). 

In the expression, (W"(ip^)H^ N \h), the presence of the function gives local- 
ization to r ~ t 1//3 , and hence, (3.3) is applicable. It follows from (3.23) and (3.3) 
that 



(3.24) 



\(W"{iP {N) )H {N \h)\ < \\H {N) \\ x i\\h\\ x i < rl/illxi +t||# (iV) || 2 Y i. 
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Combining (3.22) and (3.24), we obtain the bound 
(3.25) |£i|<* _ 1*i+*ll# W ll*i. 



Next, we compute the bound on E%. Working with (3.10), we obtain 

(3.26) \Zi\<t~ X Wx^- 

Next, we compute the Poisson bracket {E,W}. However, since {E, E} = and 
{E, M} = 0, we have {E,W} = -%{E,V q }. Substituting §J7\ and gig and 
performing numerous integrations by parts, we obtain 

{E,VM) =Rej {-i)E'^)V^)dx 

Stt R(« / (-dlij - 2r- x d r i) - 2\i)fi))(2e\d4 + aip) r 2 dr 



r=0 



?tt / |5 r V»r (d r (efr 2 ) - A9fr + ar 2 ) dr 

lr=0 S — 



v 

t -l/3 r 2 ! 



jt l/3 



8tt 



57T 



-d r (6y) + 2ar 2 ) dr 



r=0 



oo -1 
r=0 ^ £ 



t-V3 r 2 1 
.2 



<9 r (ar ) + d r (ra)) dr. 



t^r 2 1 



*l/3 



From this, it is apparent that the quadratic (in h) part of {E, V q }(ip) satisfies 
\v\- 2 \\{E,V q }&) - {E,V q }(^) - ({E,V q }'(^),h)\ 

<t 1 / 3 ||V/,||! 2 +t 1 /3||^)|| ioo || / ,||2 2+r l/3 | 



2 

L 2 



This furnishes a bound on the main term in (3.9), while E\ and £2 have been bounded 
above in (3.25) and (3.26). Thus, the bound (3.16) follows from Lemma 3.1 □ 



We next write the remainder h as the push-forward of a function /: 
(3.27) h(r, t) = e im e lvit)r/2 X(t) f{\(t)(r - q(t)),t) 

and recalling that p = A(r — q), define fi as the truncation 



(3.28) 



MM) = 0i 



Note that due to the support of 9i (which restricts \p\ <C t we have that 

A 1/2 (p/A + g)~l. 
Let 

£o = i<p, £1 = -dpV 3 , 6 = V 5 + pd p ip , £ 3 = ipcp, 
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where <p(p) = sechp, and, recalling that J = — ^i and J~ l = 2i, define 

(3-29) Kj = f (fi, |j _1 0)i 2 (-oo<p<oo)- 

Define the projection operator 

(3.30) Pfi = -K 2 {h) £ + «o(/l) 6 - «l(/l) £ 3 + K 3 (/l) &■ 

The following proposition establishes the coercivity of the functional G(h). 

Proposition 3.3. Let G(h) be defined as in (3.8), and suppose that \\h\\xi < t 2 ! 3 . 
Then 

(3.31) \\h\\ x i<G(h) + Y\ Ki \ 2 , 



i=0 



where Kj, j = 0, 1, 2, 3, are defined in (3.29). This bound is valid for e < t < to(N). 
The implicit constant is independent of e and N. 



Proof. From (3.8), 

(3.32) G(h) = (W"(^ N) )h,h)+£ 3 , 

where £3 takes the form 

£ 3 = o(t^ 3 ) J ^h 3 + 0(t A ' 3 ) J h\ 



Hence, 



(3.33) 



141 <i 4/3 ||V> (A °IU~||/i- W/bIU-II"!!^ 

< ll^ ) IU 1 ||^ll^ + t 4 / 3 ||v/ i ||i 2 ||/.|| L2 

< +rW\\h\\* xl 
1 
4 



< iS\\h\\h, 



where in the second line we have applied (3.3) in the first term and used the 3D 



Gagliardo-Nirenberg inequality (since there is no spatial localization to r ~ t 1//3 in 
this term). From (3.7), we have that W"(ip( N ^) is the operator 



1 



2 (If ^v 2 \~ 2 ) +2iv\- 2 9\d r 



W"(ipW) 

V ^ 

+ 2\~ 2 (-d 2 - 2r~ 1 d r - 4|^ {7V) | 2 - 2(^ (7V) ) 2 C), 
where C is the operator of complex conjugation. For convenience, let 

+ 2iv\- 2 e\d r 

+ 2\~ 2 (-d 2 - A\^\ 2 - 2{^) 2 C) 
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and 



D = 2[l + -v 2 \- 2 ) -2A" 2 e 



1 



Then using that d 2 r = rd 2 + 2d r , we obtain 

(3.34) (W"(i>W)h, h) L 2 m = (B(rh), rh)^^. 

Let 02 = (1 — 2 ) 1 / 2 so that /i = # 2 /i + Substituting this decomposition into the 

right side of (3.34), we get 

(3.35) 

(W"(ij {N) )h,h) LHm = (BK/j),^^)^^) + (D(r9 2 h),r9 2 h) L 2 {0<r<oo) +S 4 , 
where 

2 

£ 4 = ^ (2X- 2 {{2(d r 9 j )d r (9 j rh),rh) + 

arises from the commutator of <9 2 and 9j. These terms are lower order, however, since 
an r-derivative on h "costs" t~ 2 ^ 3 in the definition of the X 1 norm, whereas d r landing 
on Oj gives only the penalty of t~ 1//3 . Specifically, 



(3.36) 



\£*\<t^h\\ 2 xl <±6\\h\\ 2 xl . 



By (3.35) and (3.36), 



(3.37) G(h) = (B(r9 1 h),r9 1 h) + (D(r9 2 h),r9 2 h) L 2( 0<r<oo ) + £ 3 + £ 4 . 

Define the operator 

^ f (l + «9 2 + 4^ 2 + 2^ 2 C). 



Substituting (3.27), (3.28), we obtain 

(B(r9 1 h),r9 1 h) L 2 {0<r<oo) = (A/i, /i) L2( _ 00<p<oo) . 



Recalling (3.30), it is a classical fact (see for example [7 f §4]) that there exists 5 > 
such that 



(3.38) 



S\\fx ~ PfiWhai) < (Mh ~ Pfi), (h - Ph))m 



As a projection operator, P satisfies P 2 = P and P* = P (adjoint with respect to 
the (■,■) inner product defined by (3.1)). Hence, 

Wfi - Pfi\\h(u) - II/iIIl2(k) ~~ II-P/iIIl2(r)- 

Also, 

3 

iiAP/ 1 || i3(K) + iipaii^) <x;i«ii- 

fc=0 
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Applying this to (3.38), we obtain 



6\\fi\\l m <(AfiJi)l Hm +J2\ K i\ 2 - 

3=0 



Substituting back (3.27), (3.28), we obtain 

3 

(3.39) *Hi/»|li»(o<r<oo) < (Bire.h), {re.h))^^) + ^ K| 2 - 

3=0 

Directly from the definition of B, using integration by parts, we get 

(3.40) t 4/ 1^(r«||! 2( o<r<oc) < (B^h), (r9ih)) L 2 i0<r<oo) + |K/i||i 2(0<r<oo) . 
The inequalities ( |3.39[ ) and (3.40) together yield 

3 

(3.41) 5|K/>||^ 1(0<r<oo) < (B^h), {rO.h)) mo<r<oo) + l^l 2 - 

3=0 

Directly from the definition of D via integration by parts, 

(3.42) 8\\re 2 hf x i {0<r<oo) < (D(r6 2 h), (r6 2 h))v {0<r<oo) . 



Summing (3.41) and (3.42), we obtain 

(3.43) 5(||r#l/i|lxi(0<r<oo) + \\ rd 2h\\xH0<r<oo)) 



< 



{B{r9 x h), (r6»i/i)) L2(0<r<oo) + (D(r6 2 h), (r9 2 h)) L 2 {0<r<oo) + ^ \ k j 

j=0 



Combining (3.37) and (3.43), we obtain 

3 

(3.44) S(\\r0ih\\ 2 XHo<r<oo) + \\r6 2 h\\ xl{0<r<OD) ) < G(h) + £ \ Kj \* - 8 3 - £ A . 

3=0 

Again, by the fact that commutators are of lower order, we have 



(3.45) 
where 
(3.46) 



<KlkMlxi(o<r<oc) + \\ r °2h\\ xHo< , 



r<oo) I 



b\\h\\ 



+ S 5 



\£ 5 \<t^\\h\\ xl . 



Combining (3.44) and (3.45), and making use of error estimates (3.33), (3.36), and 
pll), we obtain (|3~3ll. □ 



We now come back to our original substitution (2.1) and the equation (2.2) for 
U (i.e., to our solution before the conditions (2.3) were enforced) and the equation 
(3.49) for U( N \ Consider a four parameter (Ai, &x, q±, t>i)-family of profiles: 

(3.47) 0W(Ai, qi, 6 U Vl , r, t) ^ e ie ^\(t) U^ N \\{t) (r - q(t)) , t) 
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with 

Q{r,t) = 9{t)+9 1 + {v{t)+v 1 )r/2, X(t) — AiA(t), q{t) = g(t) + q x . 

Here, the parameters Ai, q\, 9\, and t>i are assumed time- independent, and recall that 
the (time-dependent) parameters 9(t), \(t),q(t),v(t) satisfy conditions (2.3). For the 
sake of brevity, we write p = X(r—q), v = v+vi, and U^ N \ which stands for U^ N '(p,t). 
Let 



(3.48) 



FW(r,t) = idt^W + + 2|0W|V (JV) - 



Using the definition (3.47), the above is equivalent to 



FWfat) = e i0 ^A(t) (iU t {N) + X 2 U^ + 2|0W| 2 (JV) - 1 9 + y + j) <P {N) 
+ i e l& ^ (XuW + v~X 2 U^ + A(A(r - g) - Xq)Uj, N A 



+ e iB(r,t)J>^ u{ N )+i ^_ 

p + Xq P + Xq 



Recalling (2.46 ), multiplying it by A and substituting into the second and third terms 
in the above, we obtain 
(3.49) 

F^ N \r,t)= e e ^A 3 (i) 



1 1 

A2 ~ X? 



+i 



U(N) 



X ( 1 1 



A V A 2 A 5 



A(p + Ag) A(p + Ag); 



V (jv) 



+ 



xv x 9 

(l + pd p ) c/w 
2 2 



p + Ag P + Ag 



where 7 = + t> 2 /4 + -ug/2 and p = A(r — g). 
Define 



T] (r,t) = d 9l 



(3.50) 



m(r, t) 



— d 6 {N) 
X(t) 910 



Xi=l,6i=qi=vi=0 

Ai=l,0i = ( j 1 ='u 1 =0 



r j2 (r,t)= d Xl <f>W 

Ai =l,0i=gi=^i=O 

/ 73 (r,t) = (2A(t)^ 1 -g(t)A(t)^ 1 )0W 



Ai=l,0i=5i=ui=O 



Solving for id t (f)^ in (3.48) and using (3.49) (and commutation of derivatives), we 
obtain 



(3.51) 



td tVo = E"(^) Vo + H^\ 
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(3.52) 

(3.53) 
(3.54) 
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id t (r] 2 ) = E"(^) V2 - 2i\ 2 r) + Oxit- 1 ), 

= E"(^) m - 2iX 2 r ]l + Oxir 1 ). 
Proposition 3.4. Suppose that H^Hx 1 < t 2 ^ . There exists o~j satisfying 



(3.55) 
so that 
(3.56) 

(3.57) 

(3.58) 

(3.59) 



- a A < 



\Kj — a 



j\ ~ 



X 1 



\d t o- \ 



xi + \\H W \\ 



XI; 



xi + \\HW\\ xl , 



<t- 4/3 \\h\\ 2 
[X^dtXaxl^r'laol+r 2 ^ 

|Aa t A-V 2 | <r 4 /Vol + * _1 l 

These bounds are valid for e < t < t (N), and the implicit constants are independent 
of e and N. 

Proof. Let 



^(r) = 6 1 ((r-q)/q)X- l/2 r- 1 
(note that /i(r) = 0(1)) and also let 

2,M) d ^ f ^e^A^Atr-g)). 



Recall rjj, j = 0, 1, 2, 3, from (3.50) and observe that 
(3.60) Ej = Vj + L 2 m (t 1 / 3 ), 

where the error term results from cp being replaced by (see (2.44)). A change of 
variables calculation gives 



\{fll J £j) L 2 (-oo<p<oo) — \{h, J ^j) 



Define 
(3.61) 



0~ j 



then by (3.60) Kj — Oj satisfies the bound (3.55). Then 

dm = h{d t h, J _1 ?7o) + \ {h, J~ l d t r) ). 



Substituting (3.71) and (3.51), we obtain 



d t a = \{J(E'\^)h + A^\h\ 2 + 2^ N )h 2 + 2\h\ 2 h) -H^ N \ J-%) 
+ Uh,J- l (JE"(^) Vo + HW)). 
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Since J* = —J and E"(ip( N ^) is self-adjoint with respect to (•>•), we have the can- 
celation (JE"(^)h, J-^o) + (h, J- 1 JE"(ij^)r]o) = 0. Estimating the remaining 
terms, we obtain (3.56). Similarly, using the estimates (3.52) - (3.54), we obtain the 
rest of g57|-g5^ □ 

Proposition 3.5. For a sufficiently large N , independent of0<e< to(N), we have 

(3.62) \\h\\ x i < t 2N/3 , e<t<t (N). 

Proof. This will be proved by a bootstrap argument invoking the estimates obtained 
in Props. 32 3J3 and 3^ Recall that 

(3.63) \\H iN) \\ x i <t^ 2N -^ 3 . 
We make the following bootstrap assumption: 

(3.64) V < t < t (N), \\h\\ x i < 2t 2N/3 . 



Assumption (3.64) will be validated provided we can show, using Propositions 



3.2, 3.3, and 3.4, that (3.64) reinforces itself - specifically that (3.62) holds as a 



consequence. 



By (3.64), (3.63) inserted into (3.16) 



Integrating, we obtain 
(3.65) 



\d t G{h)\<t AN ' 3 -\ 



G(h) = G(h(t)) - G(h(e)) < -t 4 ^ 3 . 



By ( |3.64[ ), (|3_63j), inserted into (|3_56j), 
Integrating over [e, t] using that <x,(e) = 0, we obtain 



(3.66) 



ko|<t (27V+2)/3 - 



By (3.64), (3.63), (3.66) inserted into (3.57) and (3.58), we obtain 



lA-^Acnl <t^ 2N -^ 3 

Integrating, we obtain 



|A^A-V 2 | < t (2iV-3)/3_ 



(3.67) 



< t (2N+l)/3 



0~2 



< ± f 2N/3 

~ N 



Inserting (3.64), (3.63), (3.67) into (3.59), we obtain 

Integrating, we obtain 
(3.68) \a 3 



< J_+2N/3 

~ N 
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By (3.55), Kj can be replaced by <x, in (3.31). By (3.65), (3.66), (3.67), (3.68) inserted 
into (3.31), 



< 1^/3 
~ N 



2 

X 1 r*j 



and hence, (3.62) holds by taking N sufficiently large. 



□ 



In order to construct the solution advertised in Theorem |1.1[ we will carry out a 
compactness argument. For this compactness argument, we need higher regularity 
(H 2 ) control and tighter localization: 



Proposition 3.6. For N taken large enough so that Prop. 3.5 holds, we have the 
further bounds 



t 4 / 3 ||V 2 



, 2N i 
< t~ ~ X 



(3.69) 
and 

(3.70) !N(*)!b (R 3)<t 1/3 

for e < t < t (N), with implicit constants independent of N and < e < t (N). 
Proof. Recall 

(3.71) t h = J{E"{^)h + %^ N) \h\ 2 + 44A^h 2 + A\h\ 2 h) - ffW, 
where 

JE'\^) = +*(A + 4|^ (7V) | 2 + 2(^ {A °) 2 C) 

and C denotes the operator of complex conjugation. We have available H], control of 
h and thus seek a "higher-order energy estimate". Applying dj (j = 1, 2, 3) to (3.71 ), 
we obtain 



(3.72) 



d t djh = JE"{^)djh + A(h) + B{h, djh) - djH^, 



where A(h) contains terms of the form ip^ ■ djipN ■ h (up to complex conjugation) and 
B(h, djh) denotes terms that contain one power of djh or djh and at least one power 
of h. Thus djh, modulo error terms, satisfies at the linear level the same equation 
satisfied by h itself. Let 



Then by (|3_5j) and (j3_72j), 

(3.73) 8$ = JE\^ N) ) + JE"(i) {N) )djh + A(h) + B(h, d,h) + H {N) - djH {N) . 
Note the expansion 

(3.74) JE'$) = JE\^ N) + djh) = JE\i) {N) ) + JE"(^ N) )djh + JF, 



3D NLS CONTRACTING SPHERE BLOW-UP 



27 



where 



def 



E'{$) - £'(V> (JV) ) - E"(^)d j h 



A^ N) \d 3 h\ 2 + 2^ N )(d j h) 2 + 2\d J h\ 2 d j h. 



Using the expansion (3.74), the equation (3.73) becomes 
(3.75) d t tp = JE'$) + D-JF, 

where 

D = A(h) + B(h, djh) + HW - djH^. 

We now use the functional G that appeared in Lemma 3A_ and Prop. |3.2 , 3.3 where 
we replace ip by ip. Specifically, in place of (3.8) we take 

G(djh) = W$) - W{ip {N) ) - QV'tyW^djh). 



A slightly modified version of Lemma 3.1 follows, in which (3.75) is applied in place 
of ( 3.4[ ), and the following identity is obtained in place of (3.9) 

d t G{d 3 h) = {E, W}{$) - {E, W}(^) - ({E, W}'(^), d 3 h) - £ 1 + S 2 , 

where 

B x = {W"{^ N) )H {N \d j h) - (W'$) - W'(iP {N) ),D- JF) 
+ (W'(^ N) ),JF) 



and £2 is the same as in (3.10), but with ip replaced by ip. The proof of Prop, 
modifies accordingly, to yield in place of (3.16) 

(3.76) \dtG{djh)\ ^r^l^H^+t^- 3 

Indeed, E\ is estimated using ||W / '(V ;( ' Ar ' > )l|L 2 i$ t 1 ^-, an d 

\\D\\ L2 + t 2/3 ||l r „ t i/ 3 V£>|| L2 < t- 2 \\h\\ X i + t-*' 3 \\h\\ X x \\djhWxx 



3.2 



By an argument following the proof of Prop. |3.3 
(3.77) 



K-l 2 



i=0 



where hi = (djh, ^J" 1 ^) and H$ is given in (3.60). The equations (3.76) and (3.77) 
can be combined to yield (3.69). 

Next, we will establish (3.70). The pseudoconformal conservation law (see Strauss 

[na p. 13]) is 

1 



(3.78) 



(x + 2itV)iP\\ 2 L2(m -2t 



2t 



4 

L 4 (R3)- 
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By energy conservation, 

(3-79) Mt* m = UWIIiw - EW>) < t" 4/3 . 

By ( [2753] ), Ety) = £(?/> (iV) ( e )) ~ 1 (independently of e), and hence, g79fr yields 



L 4 (R 3 ) 



Applying this to (3.78), we obtain 



{x + 2itV)^\\ L 2 



< 



(x + 2itV)ip 



By the localization of ip( N ^ to \x\ ~ t 1 / 3 , we have 



2itV)ip 



(TV) 



L 2 



lt=e 



4/3 



L 2 



from which (3.70) follows. 



□ 



4. Proof of Theorem 11.11 



Fix iV sufficiently large so that Propositions 3J3 and 3J3 are applicable. Let to — 
t (N). For any sequence €j — y 0, let ipj be the solution to ( |1.1[ ) as described at the 
beginning of ^plwith e = 6j. By Proposition 3.5 and 3.6, we have 



ll^(^o)IU2(R3 ) + r 1/3 ||x^(to) 



\L 2 



< 1 



with implicit constant independent of j. By the Rellich compactness theorem, we can 
pass to a subsequence such that ipj{t ) — > tp in X^IR 3 ). Let ip be the solution to 
( |1 . 1 P with i/j(t ) = ip . By Proposition 3.5 



(4.1) 



sup 

6j<t<t0 



i{t)-1> 



(N) 



X 1 



< t 2N '\ 



By continuous dependence on initial conditions in the Cauchy problem, for each fixed 
< t < t , we have ipj(t) — >■ ip(t) in if 1 (IR 3 ) as j — > oo. Hence, we can send j — > oo 
in (4.1 ) to obtain 



sup ||^(t)-^ (Ar) 

0<t<i 



X' 1 



< t 2N '\ 



This yields the bound on the first two terms in (1.9). To deduce the bound on the 
third term in (1.9), we follow the argument in the proof of Proposition 3.6 utilizing 
the pseudoconformal conservation law. 
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